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£S| ' Abstract 

' The peculiar isospin properties of the b — > ccs current lead to a rich set of isospin relations 

for the B — ► D^D^K decays which are presented here. Recent high quality experimental 
data on the complete set of these decays (22 measurements) are analysed in this context, the 
isospin relations are tested and the results for the isospin amplitudes are discussed. Large 
values of the strong phases are suggested by the data. The comparison between the measured 
and expected branching fractions yields a new measurement of the ratio of branching fractions 

' Bl ?^, 4 ' S 1^~ B n - n ? = 0.86 ± 0.13. We finally discuss the implications of our findings for the 

Br(T(4:S)^B B ) ° 

measurement of sin(2/3) and cos(2/3) using these decays. 
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O 1 Introduction 

Given the difficulties in computing in a reliable and model-independent way the B meson decay 
amplitudes to hadronic final states, isospin relations are a very general and useful tool to establish 
^p" 1, relations between various B decay modes. The peculiar isospin properties of the b — ► ccs current 

are known since a long time and they have already been been used in the context of B meson 
decays The possibility that a large fraction of b —* ccs decays hadronize as B — > D^D^K 
was first suggested in Ref. 3 in the context of the discrepancy between the measured B semi- 
leptonic rate and the theoretical prediction. The same article suggested the use of isospin 
relations for the study of these decays. An additional motivation for an in-depth study of these 
channels is the possibility originally discussed in Ref. [UE1E], to measure sin(2/3) and cos(2/3) 
using these decays. Indeed they proceed through the same quark current than the gold-plated 
mode B° — > J/^K and are not CKM-suppressed to the difference of the B° — > £)(*)£)(*) modes. 

This Letter presents the complete set of isospin relations for B — ► D^D^K decays; they 
are compared to the measurements through a fit of the experimental data which determines the 
isospin amplitudes. These decays have been the object of recent experimental investigations [3 
|HJ. The last study by the BABAR Collaboration presents a complete set of measurements (22 
branching fractions have been measured) with good accuracy which is the experimental basis 
of this paper. An additional experimental complication is due to the fact that the branching 
ratios Br(T(4S) — > B + B~) and Br(T{AS) — * B°B°), needed to compare the neutral to charged 
B meson decays measured at an e + e~ machine operating at the T(4S) resonance, is not well 
known. This issue is adressed in this Letter. 
The aim of this study is manyfold: 

• verify the isospin relations using a new large set of experimental results; 
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offer some insight into the B — ► D^D^K decay mechanism from the inspection of the 
isospin amplitudes; 

Br(T(4S)->B+B-) . 



present a new measurement of the ratio of branching fractions 



Br{T{iS)^B°B a ) ' 



• discuss the implications of our findings for the measurement of sin(2/3) and cos(2/3) using 
these decays. 

2 Isospin relations for B —> D^D^K decays 

The decays considered here are B — > D^D^K, where B is either a B° or £> + , and if is 
either a if or K + . These decays proceed through a b — > ccs current through the diagrams of 
Fig. ^ Depending on the final state, the external W-emission diagram, the internal W-emission 
diagram (which is color-suppressed), or both contribute to the transition amplitude. A penguin 
diagram, shown in Fig. |2] (left plot), can also contribute to the b — ► tics current. It is expected to 
be suppressed relative to the tree diagrams of Fig. ^ and does not modify the isospin relations. 
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Figure 1: Left: internal ^-emission diagram for the decays B — > D^D^K. Right: external 
VK-emission diagram for the decays B -» D^D^K. 

The decays B° -* D^°D^°K° and B+ -* D^°D^°K + could also proceed through a 
different diagram, shown in Fig. El (right plot), which could introduce a AI = 1 amplitude. 
However this diagram proceeds through two suppressed weak vertices b — > uW and W — > su 
and a cc pair must be extracted from the vacuum, instead of a light quark pair as in the CKM 
allowed diagrams. This amplitude is therefore suppressed by at least a factor A 2 , where A is 
the expansion parameter of the Wolfenstein parametrisation. For these reasons we expect that 
AI = holds to an excellent precision. 

As already mentioned, the isospin properties of the b — * tics current are well known and 
follow from the fact that only isoscalar quarks are involved. Therefore this is a AI = weak 
transition and the final state is an isospin eigenstate. The most general expression of these 
properties is given by the relation : 

T(B+ -+ / (tics)) = T(B° - / Xtics)), (1) 

where f(tics) is obtained from the state f(tics) through a 180° isospin rotation. 
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Figure 2: Left: QCD penguin diagram for the decays B -» D^D^K. Right: CKM suppressed 
diagram with AI = 1 amplitude. 

While this relation applies to the B — ► D^D^K decays, the full structure of isospin 
relations can be obtained using the method described in [5] and summarized here. Let us 
consider a N-particle state with individual isospin quantum numbers Ik, rnt for k = 1, N (N = 3 
in our case). The isospin wave function ij){I,M) for a state of definite total isospin (I = X^k) 
can be written as 

ii(I,M)=Y,xt(I)<t>t(I,M) (2) 
t 

where t labels the invariant isospin quantum numbers (i 2 , t^..., tjv-i) of the operators T? defined 
by 

T 2 = Ii + I 2 ,... 
Tfc = Tfc_i + Ifc, ... 
I = T N ^+I N (3) 

for 3 < k < N — 1. The coefficients xt(I) do not depend on any isospin z component and the 
basis functions M) are simultaneous eigenfunctions of I 2 , 1 2 and all the T 2 .. The amplitude 
for finding the state labeled by m = (mi,wi2...mj\r) is 

<m\iP(I,M)>=J2xt(I)Umt(I,M) (4) 

t 

where 

U mt (I,M) = (h,mi;l2,m 2 \t2,mi + m 2 ) 

x (i 2 , mi + m 2 ; / 3 ,m3|£ 3 ,mi + m 2 + m 3 ) 

x ...(tiv_i,mi + m 2 ... + mN-\;lN-,m N \I, M) (5) 

and the terms on the right-hand side are Clebsch-Gordan coefficients. 

In our case, just one operator T 2 is introduced, with associated quantum numbers i 2 = 0, 1. 
The equations 0] and [3 generate the following set of relations 

A(B° - D-D°K+) = -L Al --±=A (6) 



3 



A(B°^D-D+K°) = ±A 1 + ^=A (7) 

A (B° _ D°D°K°) = -yJ\Ai, (8) 

where A\ {Aq) is the amplitude to produce the system DK with isospin quantum number 
ti = 1(0). The Ai amplitudes in these formulae are equivalent to the xt(I) coefficients of Eq. @J 
they are reduced matrix elements, in the terms of the Wigner-Eckart theorem, of the isoscalar 
Hamiltonian. 

A similar set of relations holds for charged B meson decays 

A(B 

A(B 
A(B^ 

where the A amplitudes are the same as for the neutral B decays. This isospin decomposition of 
the B D^D^K amplitudes (Eq. El to EJ) has already been presented in Ref. ^U] where it 
has been discussed in the context of tests of factorization. Identical equations hold for the other 
set of decays, B — ► DD*K, B — > D*DK and B — > D*D*K, with different amplitudes A in each 
case. In the following we have used the superscripts LL, L*, *L and ** for the B — > DDK, 
B — * DD*K, B — > D*DK and B —* D*D*K decays respectively. Equivalent relations can 
be obtained considering the isospin quantum numbers of different subsytem of the final state 
(DD, DK). The DK subsytem has been chosen here because in this case the transitions of 
Equations |H1 and ^2 proceeding only through the color-suppressed diagrams of Fig. U (left plot), 
are associated only to the A± amplitude. 

The relations presented above can be cast in the form of a triangle relation between the 
amplitudes: 

— A(B° — > D~ D°K + ) = A(B° -> D~D + K°)+A(B° -^D°D°K°) (12) 
-A{B + -> D°D + K°) = A{B + -» D°D°K + ) + A(B + -> D~D + K+) (13) 

which are depicted in Fig. 01 The two triangles for B° and B + decays are identical according to 
the isospin relations, however experimentally it is advantageous to build the triangles separately 
with the B° and B + amplitudes. 



* D°D + K°) 
- D°D°K + ) 
D-D + K + ) 



Ai, 



(9) 
(10) 

(11) 




Figure 3: Isospin triangles for the B° (left) and B + (right) amplitudes. 

We finally notice that Eq. El to El are valid not only for the total decay amplitude but also 
for each helicity amplitude separately as well as for the amplitude as a function of the Dalitz 
plot coordinates. 
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3 Study of experimental results 



The branching fractions for the charged and neutral B meson decay can be written 



Br(B+^f+) = r + /o _ N / ooi| , 3 ( / dm 2 DB dm 2 DK ) \A(B+ -+ /+)| 2 (14) 
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(2vr) 3 "32M| 

J3r(B°-/°) = ro ^^ 1 ^ (/ dm 2 DD dm 2 DK ) \A(B° -> /°)| 2 , (15) 

where r+ = 2.543 x 10 12 GeV" 1 and r = 2.343 x 10 12 GeV" 1 P] are the lifetimes of the B+ 
and B° mesons, Mb is the mass of the B meson averaged over B° and B + , m DD . and mpx are 
the invariant masses of the DD and DK subsystem, and the integral is computed numerically 
over the allowed region of the three-body phase space. In computing these integrals the small 
mass differences between neutral and charged states for the B, D* , D and K mesons have been 
neglected. 

The BABAR collaboration has recently studied the full set of B — ► D^D^K decays and 
has provided a measurement, reported in Tabled for all these modes [S]. These data, the 
most precise to date, have been obtained at the PEP-II accelerator from the reaction e + e~ — > 
T(4S) — > BB. To compute the branching fractions it has been assumed that Br(T{AS) — > 
B + B~) = Br(T(4S) — ► B°B ) = 0.5. However these equalities do not necessarily hold. In 
order to account for this factor, we have rewritten equations El and E3 m term of the rescaled 
amplitudes A = -4= where b = Br{T(AS) -» B°B°). The expression for Br(B + — > f + ) is 

then multiplied by the additional factor / + / = ^^[ri^sy^B if ) ' 

The experimental data have been fitted simultaneously using the x 2 method where the 
fitted parameters are / + /o and for each set of decays \A\\, \Aq\ and 5 = arg{A\Ao). The total 
number of fitted parameters is 13. The results of the fit are reported in Tables Q and |2 The 
overall agreement between the measured and predicted branching fractions is good as can be 
judged from Tabled Fig. Eland from the value x 2 = 8.8 for 9 degrees of freedom (ndof)- For 
this fit the statistical and systematical errors from Ref. [H] have been combined quadratically. 
This neglects the correlation between the systematical errors (common efficiencies, submode 
branching fractions, etc.). For some B decays only the sum of the branching fraction with the 
charge conjugate final state has been measured. We present in Table the fitted values for the 
individual branching fractions. 

An alternative way of displaying the experimental results and the fit results is given by the 
isospin triangles introduced above. For ease of comparison, we have normalized the triangles to 
the size of the basis (\A(B° -» D^~D^°K + )\ and \A(B + -» D^°D^+K°)\): therefore the 
lower side extends in each case from (0,0) to (1,0) and the shapes of the triangles can be directly 
compared. Given that we have only a measurement of the sides, there is a fourfold ambiguity 
on the vertex of the triangle. We have consistently chosen the same solution for its orientation. 
The seven measured triangles defined in this way are shown in Fig. [3 together with the fit result. 
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Figure 4: Results of the x 2 fit to the experimental branching fractions. The fitted branching 
fractions are shown by the stars while the points with error bars show the measured values. 
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Table 1: Branching fractions (BF) for each B — > D^D^K mode. The first error on each 
branching fraction is the statistical uncertainty and the second is the systematic uncertainty 
from The last column presents the result of the x 2 fit- 



B decay mode 


BF exp. (%) 


BF fit (%) 


B° decays through external VF-emission amplitudes 


B"- 


-> D-D°K+ 


0.17 ±0.03 ±0.03 


0.174 


B°- 


-> D-D*°K + 


0.46 ± 0.07 ± 0.07 


0.495 


B°- 


-> D*-D°K+ 


0.31±g;g| ± 0.04 


0.321 


B°- 


-> D*-D*°K + 


1.18 ±0.10 ±0.17 


1.065 


£° 


decays through external+internal VF-emission amplitudes 


B°- 


-> D-D+K 


0.08t£$? ± 0.03 


0.161 


B°- 


-> D*-D + K° + D-D* + K° 


0.65 ±0.12 ±0.10 


0.676 


B°- 


-> D*-D* + K° 


0.88±g;J| ± 0.13 


0.707 


B^ decays through internal VF-emission amplitudes 


B°- 


_> D o D o K o 


0.08 ± 0.04 ± 0.02 


0.029 


B°- 


-> D°D*°K + D*°D°K° 


0.17ig;i|±0.07 


0.181 


B°- 




0.33±gJJ ± 0.14 


0.105 


B + decays through external W-emission amplitudes 


B+- 




0.18 ±0.07 ±0.04 


0.163 


B+ - 


-» D*°D + K° 


0.41 j^f ± 0.08 


0.300 


B+ - 


_+ d°D*+K° 


0.52ig;o9 ± 0.07 


0.462 


B+- 


-»■ D*°D* + K° 


°- 78 -o:li ±0 - 14 


0.995 


B+ 


decays through external+internal VF-emission amplitudes 


B+- 


-+ 


0.19 ± 0.03 ± 0.03 


0.150 


B+ - 


-+ D*°D°K + 


0.18±g;g2 ± 0-04 


0.172 


B+ - 


-»■ D°D*°K + 


0.47 ±0.07 ±0.07 


0.459 


B+ - 


-+ D*°D*°K + 


0.53±g;lo ± 0.12 


0.660 


B + decays through internal VF-emission amplitudes 


B+- 


-» D-D + K + 


0.00 ± 0.03 ± 0.01 


0.027 


B+ - 


-+ D-D*+K + 


0.02 ±0.02 ±0.01 


0.020 


B+ - 


-» D*~D+K + 


0.15 ±0.03 ±0.02 


0.149 


B+- 


-» D*~D* + K + 


0.09 ± 0.04 ± 0.02 


0.098 
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Table 2: Results of the \ 2 fit to the experimental branching fractions. The superscripts LL, 
L*, *L and ** are for the B -► DDK, B -► 1X0*^, 5 -► D*DK and B -► D*D*K decays 
respectively. The amplitude values are in units of 10 -5 while the phases 5 are in degrees. The 
last column presents the results of the fit introducing a constraint related to other measurements 
of /+/o- 



parameter 


value 


value 


\Aj; L \ 


0.28 ±0.13 


0.25 ±0.13 


\A LL \ 


0.75 ± 0.07 


0.73 ± 0.06 


§ LL 


95 ±22 


100 ± 23 




0.27 ±0.15 


0.25 ±0.11 


\Ajf\ 


1.51 ±0.11 


1.45 ± 0.09 


5 L * 


91 ±34 


98 ±36 


\A\ L \ 


0.75 ±0.10 


0.69 ± 0.08 




1.00 ±0.11 


0.99 ±0.10 


§ *L 


111 ± 17 


116 ±14 


1^4** | 


0.71 ±0.17 


0.66 ±0.14 


1 A** 1 
\ A \ 


2.38 ±0.17 


2.27 ±0.14 


6** 


127 ± 26 


133 ± 22 


/+/0 


0.86 ±0.13 


1.02 ±0.05 


X 2 /n do f 


8.8/9 


10.4/10 


Prob(x 2 ,n do f) 


0.456 


0.406 



Table 3: Fitted values of the branching fractions for the B — > DD*K and B — > D*DK decays 
which have not been measured individually. 



B decay mode 


BF fit (%) 


B u d*-d+k u 
B° -► D~D* + K° 
B o^D* D°K 
B o ^ d°D*°K° 


0.185 
0.491 
0.160 
0.021 
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Figure 5: Isospin triangles for the B — > D^*'D^*'K amplitudes. Each panel presents the mea- 
sured vertex of the triangle, where the basis has been normalized to unity. The dotted contour 
shows the one standard deviation region. The star shows the result of the fit. We notice that 
only one triangle degenerates into a segment while in all the other cases the shape of the triangle 
presents large angles. 
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4 Discussion 



4.1 The value of f + / and the validity of isospin relations 

The value of /+/o returned by the fit is 

/+/„ = 0.86 ± 0.13. (16) 

This value is in agreement with the theoretical predictions for /+/o which lie in the 1.05-1.18 
interval |T2], as well as with other determinations of this quantity: /+/o = 1-04 ± 0.07 ± 0.04 [2] 
and /+/o = 1.10±0.06±0.05 ^3] derived from similar isospin relations for the branching fractions 
of B decays to charmonium final states. Combining these measurements obtained in B — > J j^K 
decays, rescaled using the value t + /t = 1.083 ± 0.017 we obtain f + / Q = 1.046 ± 0.056. We 
have added this constraint to the fit to the data obtaining the result shown in the last column 
of Table |2j We notice that the measurement presented here does not improve substantially the 
uncertainty on / + / and that the values and uncertainties on the amplitudes and phases do not 
change significantly using this constraint. 

The point can be investigated further. The inspection of Fig. 0] and Tabled shows that the 
branching fractions for B — > D*D*K decays deviate from the fitted values in a correlated way. 
We have repeated the fit separately for the three groups of decays final states obtaining the 
values for /+/o shown in Tabled We notice that the value measured in B — > D*D*K decays 
deviates from the experimental value in B — > J/*$>K decays by 2.95 standard deviations. 

Table 4: Values of / + /q for the different groups of decay final states. 



final states 


f+/0 


B -> DDK 

B -> D*DK + B -> DD*K 
B -> D*D*K 


1.24 ±0.43 
1.01 ±0.21 
0.55 ±0.16 



This discrepancy can be explained either by an additional systematical effect in these mea- 
surements or by a violation of the isospin symmetry for these final states. Clearly more data are 
needed to clarify this point. A high precision test of the isospin relations will only be possible 
when /+/o will be measured using a different experimental method. The large data sample 
accumulated by the BABAR and BELLE experiments will allow this measurement in the near 
future. 

4.2 Dynamical features of the amplitudes 

The amplitudes and phases extracted from the data present some distintive features. First, 
within each set, the amplitude related to the color-suppressed decays is much smaller, as 
expected. The ratios Aq/Ai are presented in Table These ratios, except for the case of 
B — ► D*DK, are close to the naive expectation |^4o|/|-^i| = N c = 3, where N c is the number of 
colors. 

Second, the central values for the relative phases 6 are in all cases close to 90°. The errors 
on these values given in Table [2] are not relevant to determine confidence intervals because of 
the non-linear relation between 5 and cos(<5) which enters the x 2 expression. To do this the x 2 
profile has been studied keeping in turn one phase 5 fixed and repeating the fit. The 90% level 
confidence intervals are 92° < 5* L < 154° and 88° < 5** < 180° while no bound can be set 
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for 5 LL and 5 L * . The superscripts LL, L*, *L and ** are for the B — > DDK, B — > DD*K, 
B — > D*DK and — > D*D*K decays respectively. From this we can conclude that there is 
a reasonable indication for large strong phases in these amplitudes. This suggests the presence 
of non-negligible Final State Interaction for these decays. This is both an important indication 
per se and has also consequences for the CP violation studies that will be discussed in the next 
section. 

Table 5: Ratios Aq/Ai from the fit to the data. 



ratio 


value 


\A^\/\A^\ 
\A$*\/\Ai*\ 
\A* L \/\A?\ 

1^0*1/1^1*1 


2.68 ± 2.44 
5.59 ±2.04 
1.33 ±0.24 
3.35 ±0.98 



4.3 Implications for a sin(2/3), cos(2/9) measurement 

All the B° —* D*D*K° are in principle good candidates for the measurement of (3. In the past 
the emphasis has been placed on the B° D*D*K° and B° D*D*K° decays (1 and 
preliminary theoretical values of the branching fractions have been presented. We notice that 
the values for the branching fractions of these modes presented in Tables ^ and El can be used 
for more precise assessments of the sensitivity of a measurement of [3 using these modes. 

In Ref.jHl, the observation of the modes B° -> D*-D* + K° and B° D~D*+K° + CC is 
reported. We notice that for B° — > D*~ D* + K° , the measured value of the branching fraction 
(0.88±g;i| ± 0.13) and the value predicted by our fit (0.707) are almost a factor two lower that 
what anticipated in Ref. [S], thereby unfortunately also reducing the comparative advantage of 
this mode with respect B — ► D*~ D* + . 

For B° -> D-D+K , Ref. jSJ reports only a 90% CL upper limit (0.17 %) which is very close 
to the fitted value 0.161 %. This means that the observation of this mode in the near future 
is possible. The estimated value of Ref. [3] (9 10 -3 ) is a factor 6 above our predicted value. 
We stress that this channel is a good candidate for CP-violation studies because of the nature 
of the final state with three pseudoscalar particles. This will facilitate the angular analysis to 
determine the helicity amplitudes. 

Finally we stress that the B° -» D*-D + K° and B° -» D-D* + K° lead to final states 
accessible by both B° and B . They can therefore be analysed in the same way as described in 
Ref. |14j . The strong phases play an important role for this analysis as the time-dependent CP- 
asymmetry amplitudes are proportional to sin(2/3 ± 5'), where 5' is the strong phase difference 
between A(B° -> D-D*+K°) and A(B° -> D-D* + K°). The possibly large values of the strong 
phases noticed above need to be taken into account for any estimate of the sensitivities of this 
analysis. 

5 Conclusion 

We have presented the complete isospin relations for the B — > D^*'D^*'K decays. These rela- 
tions have been compared to the recent experimental measurements through a fit of the isospin 
amplitudes. The overall agreement between the measured and the expected branching fractions 
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is good with the exception of a possible discrepancy for the B — ► D*D*K decays. The isospin 
amplitudes present several peculiar features which point to a dynamical origin. Large values of 
the strong phases are suggested by the data. We have also presented a new measurement of 
B B^T{Asy^B° §° ) = ^ * n a S reemen t with other determinations of this quantity. The 
implications of these results for the measurement of sin(2/3) and cos(2/3) using these decays have 
been discussed. 

6 Acknowledgments 

The author wishes to warmly thank J. Charles, J. P. Lees and L. Oliver for reading the manuscript 
and making useful suggestions. 

References 

[1] H.J. Lipkin and A. I. Sanda, Phys. Lett. B 201, 541 (1988). 

[2] CLEO Collaboration, J.P. Alexander et ai, Phys. Rev. Lett. 86, 2737 (2001). 

[3] G. Buchalla, I. Dunietz and H. Yamamoto, Phys. Lett. B 364, 185 (1995). 

[4] J. Charles, A. Le Yaouanc, L. Oliver, O. Pene, and J.C. Raynal, Phys. Lett. B 425, 375 

(1998) [Erratum-ibid. B 433, 441 (1998)]. 

[5] P. Colangelo, F. De Fazio, G. Nardulli, N. Paver and Riazuddin, Phys. Rev. D 60, 033002 

(1999) . 

[6] T.E. Browder, A. Datta, P.J. O'Donnell and S. Pakvasa, Phys. Rev. D 61, 054009 (2000). 

[7] CLEO Collaboration, CLEO CONF 97-26, EPS97 337 (1997). 

Aleph Collaboration, R. Barate et al, Eur. Phys. Jour. C 4, 387 (1998). 
Belle Collaboration, K. Abe et a/., BELLE-CONF-0104 (2001). 

[8] BABAR Collaboration, B. Aubert et ai, Phys. Rev. D 68, 092001 (2003). 

[9] M. Peshkin and J.L. Rosner, Nucl. Phys. B 122, 144 (1977). 

[10] C.W. Bauer, B. Grinstein, D. Pirjol and I.W. Stewart, Phys. Rev. D 67, 014010 (2003). 

[11] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 

[12] D. Atwood and W. Marciano, Phys. Rev. D 41, 1736 (1990). 
N. Byers and E. Eichten, Phys. Rev. D 42, 3885 (1990). 
P. Lepage, Phys. Rev. D 42, 3251 (1990). 

[13] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 65, 032001 (2002). 

[14] R. Aleksan et al, Nucl. Phys. B 361, 141 (1991). 



12 



